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Abstract 

Consider a system of particles evolving as independent and iden- 
tically distributed (i.i.d.) random walks. Initial fluctuations in the 
particle density get translated over time with velocity v, the common 
mean velocity of the random walks. Consider a box centered around 
an observer who starts at the origin and moves with constant velocity 
v. To observe interesting fluctuations beyond the translation of initial 
density fluctuations, we measure the net flux of particles over time 
into this moving box. We call this the "box-current" process. 

We generalize this current process to a distribution valued process. 
Scaling time by n and space by y/n gives current fluctuations of order 
n rf//4 where d is the space dimension. The scaling limit of the normal- 
ized current process is a distribution valued Gaussian process with 
given covariance. The limiting current process is equal in distribution 
to the solution of a given stochastic partial differential equation which 
is related to the generalized Ornstein-Uhlenbeck process. 



Key words. Independent random walks, hydrodynamic limit, current 
fluctuations, distribution valued process, generalized Ornstein-Uhlenbeck pro- 
cess. 
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1 Introduction 



In a system of independent random walks in dimension d > 1, the hydro- 
dynamic limit of particle density, under Euler scaling, satisfies the scalar 
conservation law 



with u(x, 0) = Uo(x) as the initial particle density and where t > 0, x G M. d 
and f(u) = uv, v being the mean drift of the random walks. What we mean 
by this hydrodynamic limit is the following. If we start with a sequence of 
initial particle configurations {rj^ (m) : m G Z d } satisfying 



where A C M. d is a bounded box, then the particle configurations at subse- 
quent times, r)t(-), satisfy 



where u(x,t) satisfies (11.11) with initial condition u(x,0) = uq(x). This hy- 
drodynamic limit picture indicates that on time and space scales of 0(n) we 
see shifts in the particle profile along characteristic lines of (II -ip given by 
ny + nvt, y G M. d . To capture diffusive fluctuations about these characteris- 
tic lines, we consider a box of radius 0(^/n) moving with the characteristic 
velocity v, and we look at the net inward flux of particles across the bound- 
ary of this moving box. The measurement of this net flux of particles gives 
us what we call the 'box-current' process. We generalize this box-current 
process to a distribution-valued process. In section I2.3[ we give a precise 
definition of these current processes and explain the connection between the 
current process and fluctuations about the hydrodynamic limit. 

The question we resolve in this paper is the size and the distribution of 
the scaling limit of current fluctuations. It turns out that the current process 
scaled by n~ d ^, where n is the scaling parameter, gives us Gaussian scaling 
limits. The limiting distribution-valued current process can also be attained 
as the solution of a stochastic partial differential equation. An alternate 
representation of the limiting current process as the solution of a stochastic 
differential equation is given in section 12.51 The limiting current process is 



d t u(x,t) + V x - f(u{x,t)) = 0, 



(1.1) 




x£nAnZ d 




x£nAnZ d 
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expressed as the sum of two independent stochastic integrals, where one of 
the integrals is a generalized Ornstein-Uhlenbeck process. 

In 1985, current fluctuations in a system of independently moving par- 
ticles was studied in [2] in connection with the asymptotic behavior of tra- 
jectories in an infinite particle system with collisions. Seppalainen, in [TUj, 
found the scaling limit of current across characteristics in the one dimensional 
system of asymmetric independent random walks. His work was extended in 
[7J where shifts in characteristics at the diffusive scale were allowed. In both 
cases current fluctuations were of order n 1//4 , where n was the scaling factor, 
and Gaussian scaling limits were obtained. In [TJ, a representation of the 
limiting current process as the sum of two stochastic integrals was given for 
the one dimensional case which agrees with the results in section 12.51 of this 
paper. While in the one dimensional system of independent random walks 
there is a connection between the position of a tagged particle under elastic 
collisions and a current process (see [2]), there is no such obvious connection 
in higher dimensions. 

Equilibrium fluctuations in particle density, in a system of independent 
Markovian particles in multiple dimensions, was studied by Anders Martin- 
L6f in 1976 [8]. Our results about current fluctuations hold in the equilibrium 
as well as non-equilibrium case, under the assumption that the initial dis- 
tribution of number of particles at each site on the lattice 7L d is i.i.d. It is 
interesting to note that in the absence of fluctuations in the initial particle 
configuration (i.e. if we start with a deterministic initial configuration), our 
limiting current process is a generalized Ornstein-Uhlenbeck process. Gen- 
eralized Ornstein-Uhlenbeck processes appear as the fluctuation limits of 
several infinite particle systems in the literature, for example see chapter 11 
of Kipnis and Landim [S] , PP , [1] and jS] • 

Fluctuation results for other asymmetric models in multiple dimensions 
can be found in the literature. We mention some of these. In multiple 
dimensions, equilibrium fluctuations for a tagged particle in asymmetric ex- 
clusion processes and in asymmetric zero range processes have been studied 
in [TH [TT] and [T3] respectively. The diffusion coefficient for the two di- 
mensional asymmetric simple exclusion process was studied in [T5]. In (12] . 
Sethuraman proves superdiffusivity of occupation-time variance in the two 
dimensional asymmetric exclusion process. 

In this paper we extend the results of [10] and [7J to multiple dimensions. 
We also generalize their results from a simple real-valued current process to 
a distribution-valued current process. Our proof parallels the proof in JTU] 
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and [7J in several places. Some significant differences arise in the treatment 
of a distribution-valued process. As is typical for process-level distributional 
limits, our proof is in two stages: first, we prove convergence of finite dimen- 
sional distributions and second, we prove tightness of the current process. We 
omit parts of the proof that are similar to parts in [7J and refer the reader 
to [7J for details. 

This paper is structured as follows. Section [2] contains a description of 
the independent random walks model, the current process and a statement 
of our results. The stochastic partial differential equation satisfied by the 
limiting current process and its solution are also mentioned here. In section 
E] we find the finite dimensional distribution limits of the process. We prove 
tightness of the process and complete the proof in section HI 

2 Model and Results 

2.1 Independent random walks model 

Let d > 1 denote the space dimension. We start with an initial configuration 
of particles on Z d denoted by {r] (x) : x G Z d }, i.e. r] (x) is the number 
of particles at site x at time 0. These particles evolve like independent, 
identically distributed (i.i.d.) continuous time random walks. Let X m j(t), 
m G Z d , j = 1, . . . ,r)o(m) denote the position at time t of the j-th random 
walk starting at site m. Let {rjt(x) : x G Z d } denote the occupation variables 
at time t i.e. t] t (x) is the number of random walks (particles) at site x, at 
time t. Clearly, 

r/ (m) 

mez d i=i 

x G Z d . 

The common jump rates of the random walks are given by the probability 
kernel {p(x) : x G Z d } where we assume 

The common transition probabilities of the random walks is given by 

p(x mJ ( s + t) = y\x md ( s ) = x) = Y, ^jrp {k) (y - x )> 

k=0 
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where 



P 




is the k— fold convolution of the kernel p(x). 
We make the following assumptions. 

Assumption 2.1. The initial occupation variables {rio{m) : m £ Z d } are 
i.i.d. random variables with exponential moments i.e. 

for some 6 > 0. Let p := Ei]o(m) and v : = Var?7o(m)- 

Assumption 2.2. Given {i]o{m) : m G % d }, the evolution of the random 
walks {X m j(t) — X m j(0) : m G = 1, . . . , r) (m)} are i.i.d. continuous 
time random walks on Z d starting at the origin, independent ofr) (-). We 
can assume without loss of generality that the common jump probability kernel 
{p(x) : x G Z d } is not supported on a hyperplane in M d . Let 



denote average velocity of the random walks. Define a dx d matrix a ; related 
to the second moments of the random walks, as 



— * 

v := 




X 





x=(xi,...,x d )ez d 



and choose matrix k such that 



KK 



Assumption 2.3. 



e 5lxl p(x) < oo 



x& d 



for some 5 > 0. 
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2.2 Some standard notation 

We define the norm of x = (x± : . . . , Xd) G R d as 



\x\ := max kCj . 

Ki<d 



For xGl, define 



r , f largest integer < x, if x > 
FJ := < (2.1) 
I smallest integer > x, if x < 0. 



Define 

[x] := ([ Xl ],...,[x d ]) 

for x — (x 1 ,x 2 , ■ ■ ■ , Xd) € M d . 
Let 



<9 a := 



dxi J \dxd, 

for a multi-index a = (aii, . . . , ad) with cej > 0. For any nonnegative integer 
N and multi-index a = (a±, . . . , ad), ctj > 0, define 

\\f\\ (N , a) := sup{(l + \x\) N \d«f(x)\}. 

x£R d 

Let <S(lR d ) denote the space of Schwartz functions. Recall that this space of 
functions is defined as 

S(R d ) := {/ E C°°(R d ) : ||/|| (jv , q) < oo ViV,a}. 

We denote the dual of the Schwartz space as <S'(IR d ). For any bounded subset 
B of S(R d ), let 

q B (f):= sup |/(0)|, /G5'(l d ). 

Then {g^} is a family of semi- norms on ^'(R ') which defines the strong 
topology on S'(R d ). Fix T > 0, let D([0, T},S'(R d )) be the space of mappings 
from [0,T] to iS'(IR d ) that are right continuous and have left-hand limits in 
the strong topology of iS>'(IR d ). 
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2.3 The current process and fluctuations about the hy- 
drodynamic limit 

Let 

B M := {x G R d : \x\ < M} 

be a box centered at the origin with radius M. Suppose this box is centered 
about an observer who starts at the origin and moves with constant velocity 
v. We define the 'box-current process' as the net inward flux of particles 
across the boundary of this moving box, over time. Thus the box-current 
at time t is simply the difference between the number of particles inside the 
box at time t and the number of particles initially inside the box. Scaling 
space by y/n and time by n, we construct a sequence of scaled box-current 
processes given by 

Ut, IbJ : = n~ d / 4 {Vnt(i + [nvt])-Vo(i)} 

Vo(m) 

= n ~ d/A E E [H^mAnt) E B MVE + [nvt]} - l{m E B M ^}] . 

mez d i=i 

In words, £ n (t, 1b m ) measures the cumulative net current of particles into 
the box B M ^r t up to time nt as the box moves with fixed velocity v. 

We generalize the above box-current process to a distribution valued pro- 
cess by taking the weighted average of the differences {rj nt {i+ [nvt]) —i] (i)}, 
where the weights are given by values of a Schwartz funcion. For <fi e <S(lR d ) 
and t G [0, T], define 

Ut, 0) := n- d " {<P (^r) Mi + [nvt]) - . (2.2) 

This is our current process. Observe that £„(•, •) G D([0,T],S'(W i )). We can 
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rewrite £ n (t, 0) as 




7?o(m) 



m&L d 3=1 



X m j{nt)-{nvt] 



^l{X mji (ni) =i+[nvt]}\ 



l m 



X m ,j(nt)-[nvt] \ _ i ( m \ 
^ ) ^\^)_ 

(2.3) 

In the last line of (12 .3ft we are summing over all random walks and adding 
their contributions to the current. This is a more tractable form of the 
current process compared to its initial definition in (I2.2p and we will use this 
definition henceforth for the current process. 

As in [10], here too there is a connection between the current process 
and what the author in [TU] termed the "second-order fluctuations" from 
the hydrodynamic limit. We explore this connection below. Let (Q, J 7 , V) 
be the underlying probability measure space for our independent random 
walk process. By our assumptions, under V, the initial occupation variables 
{i] (m),m G Z d } are i.i.d. random variables and are independent of the 
evolution of the random walks {X mj (-) — X m j(0) : m G Z d , j = 1, . . . f]o(m)}. 
Under these assumptions we get 



n 



-d/2 



I ~j= \ Vnt{i+[nvt}) 



-> Po 



(f){x)dx in L 2 (V) for all t > 0, 



(2.4) 

where E[r] (m)] = p . (For general initial conditions with E[t]q(x)} = Po(^t=) ; 
the right hand side would be J p(x, t)(p(x)dx where p(x, t) is the solution of a 
diffusion equation with initial condition po(-). Since po(x) = po is a constant 
in our problem, we get p(x,t) = p .) 

We can express £ n (t, 0) as the difference of two terms: 



UtA) = n d/i 



n 



-d/2 



E 



r) nt (i + [nvt]) - po <p(x)dx 
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Using the observer as a frame of reference, the first term is the fluctuation 
from the hydrodynamic limit at time t and the second term is the fluctua- 
tion from the hydrodynamic limit at time 0. Thus, £ n (i, 0) is obtained by 
subtracting the initial fluctuations from the hydrodynamic limit from the 
fluctuations at time t from the hydrodynamic limit, as seen by the observer, 
and scaling this difference by n d ^. The current process therefore looks at 
fluctuations beyond the rigid translation of fluctuations in the initial config- 
uration along the characteristic line y = x + vt. 



2.4 Results 

Recall that {p(x),x G Z d } denote the common jump probabilities of the 
random walks. Let X(t) be a continuous time random walk starting at the 
origin, with jump probabilities p(x). Henceforth, in this paper, X(t) will 
always denote such a random walk. By the martingale central limit theorem 
(refer pages 339-340 in [3]), we have 

X(nt) ~ [mTt] nB(t) 
y/n 

where B(t) is standard d-dimensional Brownian motion. The probability 
density function of nB(t) is given by 

1 v- 1 -1 

g 2t 2-<i,i a i,j XiX i 

p t (x)dx := P(nB(t) G dx) = — dx. (2.5) 

(2nt) 1 y^det(a) 

(Note: Since p(x) is not supported on a hyperplane in M. d , a is positive 
definite.) 
Define 



<7i ((s,0), (t,V)) : = / / (j)(y)ip(z)p\ t - s \(z - y)dzdy 

(2.6) 

(j)(y)ip(z)pt +s (z - y)dzdy 
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and 

a 2 (0,0), (t,V>)) := / / 0(y)*p(z)p t+s (z - y)dzdy 



<P{y)^{z)p t {z - y)dzdy - / / (j)(y)ip(z)p s (z - y)dzdy (2.7) 
+ / 4>{x)ip(x)dx 



where pt{x) is the probability density function of K,B(t) given in (12. 5p . 

THEOREM 2.1. Under Assumptions[M-{EM as n oo, £ n (-, •)-► £(v) 
in distribution in D([0,T], <S'(R )) where £(•,•) is a mean zero, distribution- 
valued Gaussian process with covariance function 

E£(s, <f>)£(t, = poo-i ((s, 0), (t, V)) + «o<t 2 ((s, 0), (t, V)) • (2.8) 



>M i 



THEOREM 2.2. Under Assumptions [U - EE3[ as n -»■ oo, ^.(-,1^ 
converges in distribution in D([0, T],R) to a mean zero Gaussian process 
with covariance as in (12. 8p where <fi = ip = 1^ 



Example. We calculate the covariance terms (I2.6P and (12.71) for a simple 
example. Take (f) = ip := 1 Bm = Yii=x 1[-m,m] and a = J dX(i . Let $„(•) denote 
the cumulative distribution function of a mean-zero real-valued Gaussian 
random variable (i.e. &,,(x) := 4 — P e~2^) and define 



i d 



111 / _ 2M 

2M{$ tt (2M) - $,(-2M)} + 2W — ( e — - 1 



Then 

O"l((s,0), = I\t-s\ ~ h+s, 

and 

o- 2 ((s, 0), (t, V)) = - h ~ Is + (2M) d . 

Remark 1. In the long-run, the current process is asymptotically a sta- 
tionary Gaussian process. This can be observed by taking the long time 
asymptotic limit of the covariance terms (12. 6p and (12. 7p . 

Remark 2. If we start with an initial configuration where ?7o(-) are i.i.d. 
Poisson(X) random variables, the system of independent random walks is in 
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equilibrium and we get 77^ (•) are i.i.d. Poisson(X) for all t > 0. Under this 
invariant distribution of particles, the limiting current process in [7] and [TD] 
are fractional Brownian motion with Hurst parameter 1/4. In the present 
paper, under this invariant distribution of particles, po = v = X and 

(2.9) 



A 



<j)(x)if>(x)dx - / / <f>(y)i/)(z)pt- 8 (z - y)dzdy 



Thus, under the invariant distribution, £(•,</>) has stationary increments, as 
E(£(t, 0) — £(s, 0)) 2 is a function of t — s. The time-indexed process £(•, 0) 
lacks the self-similarity property of one-paramter fBM. However, the process 
£(t,<f>) satisfies the following self-similarity property: 

i(at,<po Va )=a d / 2 -iit^), fora>0, 

where t] a (x) = G M. d ; this is evident from the form of the covariance 

terms ([H]) and (1277]) . 

The proof of both theorems involves two stages: first, showing conver- 
gence of finite dimensional distributions and second, proving tightness of the 
sequence of processes. We use Mitoma's theorem [U] to prove tightness of 
{£„(•,•)} C D([0,T],S'(R d )). According to Mitoma's theorem, it is suffi- 
cient to prove tightness for the sequence {£n(-,0)} C D([0,T],R) for each 
G iS(R d ). We use the tightness criteria in [2] to prove tightness of the 
sequence {£ n ( - ,0)} C D([0, T],R). A description of Schwartz space <S(IR d ), 
it's dual space iS'(!R d ) and Mitoma's theorem can be found in [5]. 



2.5 Stochastic integral representation 

The limiting current process in Theorem 12. II has the same distribution as the 
solution of the stochastic partial differential equation given below in f )2.13p . 
Define operator A on iS(R d ) (the infinitesimal generator of nB(t)) as 

1 d 

M(x) ■= g Yl *iAi<t>i x )- ( 2 - 10 ) 
Let {T t ,t > 0} denote the semigroup associated to A. Then 

T t (4>)(x) = / (t>(y)p t {y - x)dy. 

JM. d 



11 



Define covariance functions 



» d 

Qi(0,V):= / y~] aijdi<f>(x)djif)(x)dx } (2.11) 
jRd i,i=i 

and 

Q 2 {(f),^):= I A4>{x)A^{x)dx. (2.12) 



Consider the following spde, 

dZ(t, •) = sffrdW; + AZ(t, -)dt + y/v^JF(-)dt] Z(0,-) = 0, (2.13) 

where Wt is a centered <S'(IR d )- Wiener process with covariance Qi(-.-), {J-"(0), G 
«S(R )} is a Gaussian random field on iS(M d ) with covariance Q 2 (-, •) and J 7 
and W t are independent of each other. We can write 



Wf= [ V0(x) • KWt{dtdx) (2.14) 

where W\(t, x) is <i-dimensional space-time white noise on M + x M d . 

= / A<j){x)W 2 {dx) (2.15) 

where W^O^) is a one- dimensional white noise on M. d independent of the white 
noise W%. Observe that the stochastic differential equation (I2.13P has two 
independent sources of randomness and thus we get the solution to ( 12.1 3p to 
be the sum of two independent stochastic integrals 



Z{t, 0) =^p- [ V (T t _ s <j)(x)) ■ KW x {dxds) 

+ I I AT t _ s (f>(x)W 2 {dx)ds. 

J0 JR d 

Denote the two stochastic integrals in ( 12. 16f) as Zi(t, 0) and Z 2 (t,(f)) so that 

Z(t, 0) = y/p^Z^t, 0) + y/v^Z 2 (t, 0). 

The first term v /p^Z i (t, 0) represents space-time noise created by fluctuations 
in the random walks. The second term y/vo~Z 2 (t, 0) indicates the fluctuations 
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in the initial configuration that get propagated as the system evolves. In the 
absence of fluctuations in the initial configuration, i.e. if we take the initial 
configuration to be deterministic (vq = 0), the second source of randomness 
in the spde (I2.13p . y/voF^dt, disappears. The spde (I2.13P in this case then 
reduces to 

dZ(t, •) = y/frdWi + AZ(t, -)dt, 

where the solution 

Z(t, <j>) = y/p^f V (T*- s 0(z)) ■ nW 1 (dxds) 

is a generalized Ornstein-Uhlenbeck process jl] . 
For t,se [0, T] and <f>,ij) e S(R d ), 

Cov(Z(t, 0), Z(s, i(>)) = po Cov(Z 1 (t, 0), Zx(s, ^))+v Cov(Z 2 (t, 0), Z 2 (s, </>))• 

Elementary calculations show that 

Cov^t, 0), Z^s, i))) = ((s, 0), (t, 

and 

Cov(Z 2 (t, 0), Z 2 (s, = a 2 ((s, 0), (t, V)) • 

We can conclude that the solution of the spde and the limiting current pro- 
cess in Theorem 12.11 are equal in distribution as they are both mean- zero, 
distribution-valued Gaussian processes with the same covariance. 

3 Convergence of finite dimensional distribu- 
tions 

We begin the proof of Theorems 12.11 and 12.21 with showing convergence of 
finite dimensional distributions. 

Note: Henceforth, in this paper, 'c' will denote constants that change from 
line to line in calculations. 

Fix iV G N and choose (U, fc) G [0, T] x S(R d ) for i = 1, . . . , N such that 
(tj, (pi) 7^ (tj, <pj) for i 7^ j. Let (#!,..., 6 N ) G WL N be an arbitrary vector. 
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LEMMA 3.1. As n — > oo, Yli=i @i£,n(ti,4>i) converges to a mean-zero Gaus- 
sian random variable with variance 

d 

o- 2 = ^2 did i (P° ai ((*<> &)> &)) + WoCT2 ((*<> &)> &))) • 
i,i=l 

Proof. Define 

??o(m) 

U m (t, 0) := n- rf / 4 £ [0 - (^)" . 

j'=i 

Let W m = Y,i=i8iUm(ti,4>i)- Denote U m (t,4>) = U m (t,4>) - EU m (t,</>) and 
= W m — EW m . Choose r(n) so that r(n) — > oo slowly enough that 

(r(n)) d • E[r/o(x) 2 l{r7 (x) > n d/8 }] -> as n ->■ oo. (3.1) 

Since {^o( m ) : m £ are i-i-d. , a simple calculation shows that Erj t (m) = 
Po, Vt > 0, m G Z d and hence E£ n (t, 0) = 0. Rewrite 

AT AT 

i=l 8=1 meZ d 

This sum can be split into the following two sums: 

N 

i=l |m|<r(n)-y/n |m|>r(n)v^i 

We can now apply the Lindeberg-Feller theorem to the sum EimkrfnJVn^ 
and show that it converges in distribution to the mean zero Gaussian random 
variable with variance indicated in Lemma I3TT1 The proof follows in the same 
vein as the proof of Lemma 2 in [7] and we skip it. 

We next show that the second sum in (I3.2p converges to in L 2 as n — > oo. 
The proof of this differs from the analogous step in [7] and so we give the 
proof in detail. In [7J we use large deviations to control contributions to the 
current from distant particles. This is not enough in the present situation and 
we need to appeal to the rapidly decreasing property of Schwartz functions 
to bound these contributions. 
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N 



By Schwarz inequality, 

/ V 

E \ E W A £H»H 2 E E EUAuAif 

\jm\>r(n)^ri J t=l |m|>r(n)\A* 

where ||0|| 2 = X^li^i- It suffices to show that for a fixed t > and 6 
S(R d ), 

\m\ >r(n)y/n 

as Tl — t- OO. 



^ £f/ m (t,0) 2 

|m|>r(n)y / n 

< cn"^ 2 £ E 

\m\>r(n)y/n 



Vo(m) 

0^ X m ,j(nt)-[ngt] - 



n 2 



^m,i("t) — [nvt] ■ 



+ ™ d/2 E * 2 ($) 



|m|>r(n)\/n 

since (a + 6) 2 < 2a 2 + 26 2 and rj Q (-) have bounded moments by Assumption 



< cn- d / 2 



E w 



2 f Xm,i(nt)-[nvt] - 
v 7 " 



+ E ^) 



|m|>r"(n)y/n 



|m|>r(n)v^n 



since f]o{-) have bounded second moments and are independent of the random 
walks 



cn 



-d/2 



E 



E U 2 ( X -'"^' [ "' , " )l{AVi(»t) 6 B„ ( „)^ + [nvt]} 



\m\>r(ri)y/n 
k 2( Xm,i_{nt)-[nvt] 



+ E U\ A - 1{n X [nmi )l{X m ^nt) i B M{n)VE + [nvt]}) + 2 (^ 



(3.3) 



where M(n) is chosen so that M{n) — > oo as n — > oo and M(n) = o(r(n)). 
= cn- d / 2 (h + I 2 + I 3 ) 
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where 7 1 ,7 2 ,/3 denote the sums of the three terms resp. in the expression 
above. We show below that the scaled limits of each of these sums goes to 
as n — ?■ oo. 

n d l 2 h < cn- d ' 2 P (X(nt) G B M(n)VE + [nvt] - m) 

\m\>r(n)y/n 

since is bounded 



= cn- d ' 2 P ( X ( nt ) ~ Ml i B r[n) ^ + i) 

<cn- d ' 2 P (X(nt)-[nvt]^B hr{n)VE ) 



as \r{n)^/n < r(n)y/n — M(n)y/n (recall that M{n) = o(r(n)) 



cn- d ' 2 



£ ^(l^(^)-M]|>^(n)v^) 

! M(n 

< cn 



^ dl2 E\X{nt)-[nvt]\> 



(r(n)y/n) 
by Markov inequality, r > 

= cn r/2 n" r/2 M(n) d (r(n))- r 
(as £|X(nt) - [ra/t]| r is 0{n r l 2 )) 
->■ 

as n — )> oo by taking r > d. 
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The following requires the rapidly decreasing property of Schwartz func- 
tions. 



^ 2 I 2 <n- d / 2 E E (< 

meB c . . _ L>M (n) 



,/ " "/, < „ " " > > t[d 2 ( X mil (nJ)-[nvtV 



meB c . . _ L>M (n) 

since G iS(M d ), \4>(x)\ < CjX~i for any j G N 



x l{X m ,i(nt) G 5^ n + M]} 

L 2 " 



<^ d/2 E E 4W-M]e%nB w - m 



< n- d / 2 EE E jij p ( x ( nt } - m = % ~ m ) 

L>Af(n) 

by summing over m first and then % 

L>M(n) L>M (n) 

as n — 7- oo by choosing 2j > d and since M(n) — > oo as n — > oo. And finally, 

" d/2 / 3 < c / 2 (x)dx -> 

J \x\>r(n) 



n 



as n — » oo by a Riemann sum argument and since r(n) — >■ oo as n — > oo. 
Thus the right hand side of (13. 3p goes to as n — >■ oo and since 



e\ J2 Wm <ii#ii 2 E E E Um{u, 

\m\>r(n)y / n J *=1 |m|>r(n)\A* 



iV 



we conclude that the second sum in (13. 2p goes to in L 2 . This proves Lemma 
I3TT1 □ 
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4 Tightness and Completion of the Proof 



Choose a > and 

/3>d/A + da + l. (4.1) 

Fix E 5(M d ). We prove 

PROPOSITION 4.1. £„(•,<£) in D([0,T],R). 

To prove Proposition 14.11 we use the tightness criteria in [2] . In sections 
14. II and [4. 21 below, we check the two tightness conditions of Proposition 5.7 in 
[2] for the sequence of processes {£«(•, 0)}. The first condition is a moment 
bound condition and the second involves the modulus of continuity. 



4.1 Verifying the first tightness condition 

Let s, t E [0,T] and without loss of generality we assume t > s below. Let £ 
denote the centered current process. Choose an integer 

r > 2(3 > 2. (4.2) 

We show that 

LEMMA 4.1. When (t - s) > n~ & , 

E \Ut, 0) - Us, <t>) T < C r (t - sY (4.3) 
where o > 1, and C r is a constant depending only on r. 
Proof. Let 

770 (m) 

A •— ^0 ( Xm ' J ( - nf - ) ~ ) ( Xm < 3 ( ns )~ [ ngs l ) ' 

We can write 

£n(t,0)-£n(s,0)=n- d / 4 ]T A 



Then 



2r 



(4.4) 



where A m = A m — EA m . 

To appropriately bound f !4.4p . we first find the following moment bound. 
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LEMMA 4.2. For any integer 1 < k < 2r, there exists a constant C that 
depends on r, such that 



E \ A A < C{n- 1/2 + )ri 



d/2 



Proof. 



\ k <2 k 



by applying (a + b) < 2 (\a\ + \b\ ) and then Holder's inequality 



by Jensen's inequality and independence of 770 (•) from the random walks 



< c2 



2r 



(4.5) 



since r)o(m) has bounded moments. 
Define 

:= . . . , <9 d 0(x)) , 

where <9, denotes the partial derivative with respect to the ith co-ordinate of 
x. Let ip(x) := (1 + l^l) - ^ for some positive integer N. 



jjj \^^ X mt i(nt)-[nvt] ^ _ ^ X m ,i(m)-[nw] 



^/,^ X m ,i(ns)-[ngs] 



^.lM-liift] _ ^m.i ("■«)- [ni/s] 



^ m jy m ,i(nt)-[ngf] I ^ I X m ,i(ns)-[nt?s] 1-1 



Xl{| 



-Xm,l(ns) — [ngs] 1 I X m; i(nf)-[ni?f] 



v 7 " 



> 



•y/n 



n 1/ 



19 



by applying the mean value theorem and \D(f>(x)\ < cn^(x). 

We will treat one of the terms in the above inequality, the other being 
similar. 



C N 

m<=Z d 



E 



Xm,i(ns) — [nvs] 



X 1 



Xm,l(nt) — [nvt] Xm,i(ns) — [nvs] 



y/n 



X mt i(nt)-[nvt} 



> 



Xm,i(ns) — [nvs] 



[l + L) kN 
x 1 



X m ,i(nt)—[nvt] X m ^(ns)—[nvs\ 



sfn 
(ns)—[nvs] 



eB c L nB {L+1) ) 



+ L) 



kN 



E 



X{nt)-X{ns)-[nvt] + [nvs\ 
s/n 



x 1 |x(ns) G B C LVE n B {L+l)VE + [nvs] - m} 



L>0 

by summing over m 



X{nt)-X{ns)-[nvt] + [nvs\ 
s/ri 



< C N 



-kN+d-l 



L>0 



n 



d/2 [{t _ s) k/2 + n -k/2 li 



using central limit theorem 



< cn d/2 [Vt^ + 



n 



-1/21 



by choosing N large enough so that the sum in parenthesis is finite and by 
bounding (t - s) k/2 < (2T) ^ ^/t^~s and n^ 2 < n^ 2 . □ 



LEMMA 4.3. There exists a constant C depending on r such that 

E ( Yl A ™) ^ c ^ + ( n_r/2 + (* - s ) r/2 ) nrd/2 > 

\m<=Z d / 



(4.6) 
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The proof of the above lemma uses Lemma 14.21 and is along the same 
lines as Lemma 8 in |7|. (In fact Lemma 8 in [Tjj is a special case where r is 
6.) We omit the proof. The reader is referred to Lemma 8 in [7] for the idea 
of the proof. 



By (I4.4p and the above lemma we get, when \t — s\ > 



n 



e \ut, 0) - us, 0) r < c i n ~ rd/2 + n ' r/2 + (* - s Y /2 ] 

< c \{n~ p )^ + (rT p )^ + (t- s) T/2 

< c(t - s) a 

where a > 1 by our choice of r > 2/3 > 2 in (14. 2 p and where c is a constant 
depending on r. □ 



4.2 Verifying the second tightness condition 

The second tightness condition involves proving that 

lim P(w^(n" ) > e) = 0, (4.7) 

where 

WiS n ~ P ) ■= SUp |^ n (*,0) -£nO,0)| 
\t-s\<n~P 

is the modulus of continuity. 

Recall the definition of the step function [•] in (12.11) . For each n, we will 
divide [0, T] into subintervals such that [nvt] is constant in each subinterval 
and each subinterval has length less than . We construct a family X n of 
left-closed right-open subintervals by following the steps below. 
Step 1: Set a = 0. 

Step 2: Let b\ be the largest real number greater than a such that [nvt] is 
constant as t varies over interval [a, bi). 
Step 3: Set b = min{&i, a + n _/3 }. 
Step 4: Include [a, b) in Z n . 

Step 5: If 6 > T then stop, otherwise set a = b and repeat from Step 2. 

Clearly X n forms a minimal covering of [0, T] that satisfies the following 
two conditions: 

1) V[a, 6) G X n , |6 — a\ < rT 13 and [nift] is constant in [a, 6), ^ 

2) the subintervals in X n are disjoint. 
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A little thought gives us 

\L n \ < d\TwF + 1] for large n 

(i.e. number of such subintervals is at most d\Tn" + 1] where d is the dimen- 
sion) . 

To verify the second tightness condition, it is sufficient to show 
LEMMA 4.4. For any < T < oo and e > 0, 

lim V P(sup |£ n (t,0)-£ n (a,0)|>e) = O. (4.9) 

The reason that Lemma 14.41 is sufficient to prove (14.71) is the following. 

For any < s < t < T with t — s < n~@, there exist [ai,&i), [02,^2) £ 
such that s G [ai,&i) and t G [a 2 ,& 2 ) ([ a i)^i) m &y be equal to [02,62))- By 
the triangle inequality we get 

|£n(t, 0) - 0)1 < 0) - £n(>2, 0)1 + lfn(a 2 , 0) - £ n (a 2 + n"", 0)| 
+ |fn(a2 + n _/3 ,0) -^ n (ai,0)| + |£„(a x ,0) -£„(s,0)| 

(4.10) 

By applying Lemma [4.11 and \X n \ < d\TvP + 1], we get 
P{w u {n^) > e) < 2 P ( SU P ^) - &>( a > $1 > e / 4 ) + C ^ (1 " CT) 

r i.\^t t£\a,b) 

[o,6)ez„ 1 y 

As n — >■ 00, the second term in the above inequality goes to as o > 1. 



Proof of Lemma \4-4\ Let 



(t):= P\ suv\UtA)-Ua,<f>)\>£)- (4.11) 

[a,6)6Z n V* 6 ) / 

We want to show that 

lim (f) = 0. 

We start by giving a general idea of the proof. The current process is written 
as the sum of two parts: contributions from particles starting inside the box 
B n i/2+ a (recall a > from beginning of this section) and contributions from 
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particles starting outside the box B n i/ 2 + a . In the first case we bound the cur- 
rent by the number of jumps executed by particles initially within the box 
B n i/2+ a . In the second case we conclude that particles starting outside the 
box B n i/2+ a can either travel a distance of order greater than y/n towards the 
origin, which would be a large deviation, or they would remain sufficiently 
far from the origin, in which case the value of <fi at that distance from the 
origin would be small. 

Choose 7 such that 

< 7 < a. (4.12) 

We first find a bound on the expected number of particles starting outside 
the box B n i/2+ a , which enter the box _B n i/2+ 7 at some point in time interval 
[0,nT\. 

LEMMA 4.5. Define 

VO (m) 

Ni := J2 J2 l{X mJ (nt) G 5 n i /2+7 + [nvt] 

| m | >n l/2+a j = l 

for some < t < T}. 
Let I > 1 . Then there exists a constant c independent of n such that EN\ < 

cn d/2-la+d'Y 

Proof. Recall that X(nt) denotes a random walk starting at the origin, with 
the same distribution as the evolution of the particles. Let X(nt) := X(nt) — 
[nvt] . 

ENi = p P{X{nt) e B n i/2+-, - m for some t G [0, T]) 

| m | >n l/2+a 

since r] is independent of the the random walks 

= p Yl p (X(nt) (£ B n i /2+a + i for some t € [0, T]) 

*e- B „i/2+ 7 

< p p ( x (nt) i Bi nl/ 2 +a for some t e [0, T]) 
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(4.13) 



for large enough n, \w}l 2JrOL < n l / 2+a — n 1 / 24 ' 7 , as 7 < 



a 



Po 



J2 p ( ™p \X(nt)\ > -n 1/2+a ) 



< cn d/2 (n f ) d ^- 

( —71 1 



0<t<T 



E\X(nT)\ l 

(l n l/2+a)J 



by Doob's inequality for any I > 1 

< cn d l 2 {n^) d n l l 2 n la - l l 2 
as E\X{nT))\ l is 0(n'/ 2 ) 

< cn d/2 - la+d \ 

Let [a, b) £ X n . Define 
We can rewrite (f) in (14. lip as 



X m j(na)—[nva] 



(t)= E P 



[a,b)eX n 



sup 

t£[a,b) 



> 



+ J2 p su p 



[a,6)eX n 



r?o(m) 

E E^'w 

m6B„i/2+ J =1 

770 (m) 

E E ft Sw 



n d/4 e 



te[o,6) 



> 



n d/4 £ 



□ 



(4.14a) 
(4.14b) 



The two sums reflect the split in contributions to the current process from 
particles starting within B n i/2+ a in (I4.14a[) . versus particles starting outside 
B nl/2+ a in fOibj) . 

We first show that the second term f l4.14bl) goes to as n — > 00. To do 
this we split ( 14.14bl) into two sums, the first containing contributions from 
particles that enter B n i/2+~, at some time in [0, nT] and the second containing 
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contributions from particles that never enter B n i/2+~, 

n d/4 e 



sup 



7?o(m) 

E EOo 

meBc n i/2 +a i =1 

Vo(m) 



E p 

[a,b)ez n \ te ^ 

^ ^ cl{X mJ (nt) G S nV2+7 + [nvt] 



[o,b)GJn Xm es- 1/2+Q j=l 



for some t G [0,T]} > 



n d/4 e 



(4.15a) 



+ E K sup E EhSw 

[a,6)eX„ V * e[a ' b) '»6^ 1/2+Q j=l 



1 {X mJ (nt) G ^ 1/2+7 + [nvt], Vt G [a, 6)} > 



n d/4 e 



(4.15b) 



To get ( 14. 15 aft , we bounded h^j by some constant c times the indicator 
function, since is a bounded function. 
Now using Lemma I4"3| we get 

d/A, 



/ n d/4 e\ 

w < Yl p ( cNi -—) 



by Markov inequality and since \I n \ < d[Tn^ + 1] 

< cn^n- d / 4 n d / 2 - la+d ^ 

by Lemma \A. 51 Choose T large enough so that + d/4 — la + d'j < 0. Then 
the right hand side — > as n — )■ oo. 

We use the property that Schwartz functions are rapidly decreasing to 
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show (14.15bl) goes to 0. 

tijo(m) 
E E 
meSC l/2+a j = l 



c n {1 + L)~ n i(l< M | W^)-M1 | < (L + l)\ > 



L>nT 

since |0(x)| < Cw(l + |^|) _Ar 



£ E E E Mi + vr" 

[a,6)eZ n m£B^ 1/2+a L>r^ 

P ( X - l( ^" [wgt] € 5£ n % +1) for some t G [a, 6) 



by Markov inequality and since r/o is independent of the random walks 

= ^"« E E E*d+£)- N 

[a,6)eXn m£B= 1/2+a L>nT 

x P ^X(nt) — [raft] G n ~~ m ^ or some * e K 

<^ Cw Yl E( 1+L ) _JV ( L ^) d_1 ^ 

[a,b)eX„ L>nT 

by summing over m and allowing cjy to absorb all the other constants 

„d/4+/3 

< C7V ^- Yl L ~ N+d ~ l 

L>n~< 

since |X n | < d[Tn^ + 1] 

< c N n d/i+p - l(N - d) 

by choosing iV large enough so that d/4 + /3 — j(N — d) < 0. 

We finally turn to ( 14. 14a|) . To prove (I4.14aj) goes to 0, we show that the 
number of particles initially inside B n i/2+ a that jump during a time interval 
of length n 1_/3 , is stochastically smaller than 0(n d//4 ). 
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Recall, by definition of X n , [nvt] is constant for t G [a, b) G X n . Therefore 
for t G [a, b), 

\h [ ^(t)\<Cl{G^h 

where G^ b - := {X m j(nt) jumps during time interval t G [a, b)} and C = 
2 sup 2 



(x). Define 
A^) :=\ogE 



rjo (m) 
m e B n l/2+a i = i 



and 



A^) := logS [exp {i{g£J } }} 



(4.16) 



Let I1(?t, 1 be a Poisson(n 1 ^) random variable. 

^M) = i g|i + (e — 1)P (X m) i(nt) jumps during time interval t G [a, 6))} 

< log {l + (e - i)p (n(n 1_/? ) > 1) } 



since b — a < n 13 

< (e-l)P(U(n 1 -^) > 1) 

by Markov inequality. 

by independence of tjq from the random walks 

< J2 log[l + p cn 1 - /3 (e-l)+o(n 1 - /3 )] 

mEB nl/2+a 

by f)4.16p and since rj (m) has exponential moments 

< cnl ~ P ^ C n d ' 2+da+1 -P 

m&B nl/2+a 



(4.17) 
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Putting all this together, we get 



wi < E p E E > > 

[a,b)£X n \m6B n i/2 +Q J=l 

since |^(t)| < C1{G [ $} for t G [a, b) 



~2C~ 



< 



cxp 



[a,6)eX„ 

by Markov inequality 



_ n d/i e 
2C 



E 



Vo(m) 

E E 



< 



E exp 

[a,6)6Z„ 
< d[Tn p + 1] exp {- 



. n d/4 e 
2C 



exp 



{A M) } 



-cn 



d/A 



(1 



<2/4+d»+l-/3 



)} 



by applying fTCTTj) and \X n \ < d[Tn p + 1] 



->■ 



as n — )• oo, since /3 > c?/4 + da + 1 by definition (14. ip . 

This proves Lemma I4~4l and thus verifies the second tightness condition ( 14. 7ft . 

□ 



Proof of Proposition \4-l\ Lemma 14.11 and ( 14 .TP satisfy the tightness criteria 
in Proposition 5.7 of [2], thus we get £ n ( - ,0) is tight in D([0, T],R). □ 



Proof of Theorem \2.1[ We invoke Theorem 4.1 in [9] which states that Propo- 
sition ETTl is sufficient to prove that the sequence {£ n ( - , •)} i s fight in D([0, T], S'(M. d )). 
This, together with Lemma 13.11 proves the theorem. □ 



Proof of Theorem {KM The proof of this theorem follows almost verbatim 
from that of Theorem 12.11 A few places where the proof differs from that of 
Theorem 12.11 are highlighted below. 
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When proving tightness, the proof of Lemma T4.2I is different. From equa- 
tion ( 14. 5 p onwards in Lemma 14.21 the proof differs as follows. 



E \ A m\ k < c2 2r 



k 

' X m ^(nt)-[nvt] \ _ , / X m ^(ns)-[nvs] 



< c2 2r l p ( c m) + P(D m )} 

where C m = {X mA (nt) e B My/E + [nvt], X mjX {ns) <£ B My/E + [nvs]} 
and D m = {X m>x {nt) B My ^ + [nvt], X mA (ns) G B M ^ + [nvs]}. 



P(Cm) = P{X{nt) e B MVE + [nvt] - m, X{ns) <£ B M ^ + [■ 



nvs] — m) 



= J2J2 P ( X ( nt ) - X ( ns ) e B m^ + M*] - [nvs\ - j | 

X(ns) = j + [nvs] — m) x P(X(ns) = j + [nvs] — m) 
= P(X(n(t-s))eB M ^+[nvt]-[nvs]- 3 ) 

= P{x n (t,s)eB MVTi + j) 

where X n (t, s) = X(n(t — s)) — [nvt] + [nvs] 

= J2P(Xn(t,s) = k) HkeB MV7l + j} 

k& d H B Myfr 
My/n 

< J2 E P(X n (t,s) = k)n^\k\+ Yl P(X n (t,s) = k)cn d / 2 

1=0 k£Z d :\k\=l \k\>My/n 

< n^ E\X n {t, s)\ + cn d/2 P(X n (t, s) > My/n) 

since E \X(n(t — s)) — n(t — s)v\ is O (^Jn{t — s) j 

< c[n^(Vn~(F^sj + 1) + n d l 2 U{t X ] < ^[v/^ + n- 1 / 2 ]. 
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Similarly, 



The proof of (14 .TP holds for indicator functions without any modifications. 

□ 
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